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ONE  CONTROL  ALGORITHM  FOR  FINAL  RATE  OF  DESCENT  OF  AUTOMATIC  VEHICLES 
IN  MARTIAN  ATMOSPHERE 


N.  H.  Ivanov,  A.  I.  Martynov 


ABSTRACT 


Proposed  is  a staple  control  algoritha  for  the  final 


rate  of  descent  of  antonatic  vehicles  in  the  Martian  atnosphere  which 
gives  a ainiaun  speed  at  a certain  final  altitude.  The  lifting  force 
vector  is  controlled  by  changing  effective  quality.  Presented  here 
are  nanerical  results  obtained  froa  estiaating  the  effectiveness  of 
the  proposed  algoritha  for  two  hypothetic*!  leading  vehicles  having 
the  sane  value  of  available  quality  ■ 0.3  but  a different  reduced 
load  on  the  frontal  surface:  p^  = 80  kgf/a*  and  Px  * 250  itgf/a*.  BHD 
ABSTRACT 
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The  rarefied  Hartian  ataosphere  coaplicates  the  lauding  of  space 
vehicles  on  its  surface  when  aerodynanic  braking  is  used.  In  a nuaber 
of  studies  (see,  for  exaaple  [ 1 }-[  3 ])  it  has  been  shown  that  only  in 
the  presence  of  autoaatic  guidance  facilities  assuring  relatively 
accurate  entry  into  the  ataosphere  (navigation  corridor  of  entry  with 
respect  to  altitude  of  conditional  pericenter  bHr.  » ISO  ka)  and  in 
using  landing  vehicles  (LV)  with  lifting  force,  can  alaost  coaplete 
daaping  of  the  energy  of  the  LV  be  achieved  by  using  the  ataosphere. 
Here  one  of  the  central  probleas  is  that  of  creating  a systea  for 
controlling  landing  (LCS).  The  considerable  scatter  in  entry 
conditions  (priaary  with  respect  to  altitude  of  the  conditional 
pericenter  of  entry  trajectory  # the  high  degree  of  indeterninacy 
in  the  value  of  ataosphere  paraaeters  [2]  in  coabination  with  the 
requireaent  for  aaxiaal  siaplicity  and  reliability  of  the  LCS 
coaplicate  the  solution  to  this  problea. 

The  problea  of  controlling  a LV  in  the  Martian  ataosphere,  at 
least  in  the  first  stage,  can  be  foraulated  as  a problea  of 
controlling  the  rate  of  descent  V*  at  e given  final  altitude  H,.  This 
rate  should  be  ainiaal  for  creating  the  aost  favorable  conditions  for 
operation  of  a soft  landing  systea  (SLS).  It  is  assuned  that  shea H=H% 
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the  refiae  of  steady-state  descent  has  not  yet  begun.  Plight  range 
L plays  no  substantial  role,  and  can  be  ignored.  Hith  automatic  LV  in 

■ind,  load  factor  liaitations  need  not  be  considered. 

4 < 

Given  in  this  article  is  an  analysis  of  an  automatic  control 
algorithm  for  final  speed,  which  has  a sisple  structure  and  requires 
only  sisple  Measuring  and  calculation  devices  for  its  achievenent.  At 
the  saae  tiae  the  algoritha  does  possess  the  necessary  universality 
and  flexibility  for  use  under  different  conditions  of  approach  to 
Ears.  The  proposed  algoritha  provides  reliable  and  sufficiently 
accurate  guidance  of  the  LV  to  the  assigned  final  altitude  H = HV 
with  a velocity  value  close  to  the  ainiaal  while  aaintaining  the 
altitude  liaitatioa  /Y>  Studied  in  this  case  was  the  equation 

using  the  angle  of  roll  at  which  the  effective  aerodynes!?  equality 

5 K^  — Kt  cos  7, 

where  Kt  — Kpttn  is  the  value  of  quality  at  a balanced  angla  of  attack, 

T - of  roll. 

Roraal  control  prograa.  studies  conducted  earlier  [3]  showed 
that  results  close  to  optiaal  froa  the  standpoint  of  obtain!*?  V* 
for  LV  with  average  values  of  reduced  load  on  the  frontal  surface  p^ 

* G/Cg S,  provides  a control  prograa  with  a single  switching  of  the 
roll  angle.  The  LV  enters  the  flartian  ataosphere  with  a ainiaal  value 


DOC  * 0196 


PAGE  4 


of  effective  geality  it  a certain  aoaent  in  tiae 

there  occurs  the  transition  to  flight  with  — Km„ (t  = 0).  thereafter 
vehicle  aoveaent  occurs  with  a constant  value  of  K*.  Despite  the 
fact  that  the  position  of  the  switching  point  is  largely  dependent  on 
aany  factors  (priaarily  on  the  paraaeters  of  ataosphere,  initial 
entry  conditions  and  characteristics  of  the  vehicle)*  as  a whole  the 
contrpl  prograa  is  extreaely  siaple  and  can  therefore  be  ased  in  the 
capacity  of  a noainal  prograa  for  creating  a control  algor itha.  For 
this  purpose  ve  aust  deteraine  by  onboard  Beans  a switching  point 
which  will  satisfy  the  conditions  of  the  problen  which  has  been  posed 
for  obtaining  V««i„  with  the  actual  paraaeters  of  the  trajectory* 

IV,  and  the  ataosphere. 

Selection  of  switching  line.  In  the  general  case  the  switching 
point  for  each  specific  i-th  trajectory  can  be  assigned  by 
coordinates  Hta,  share  • is  the  angle  of  inclination  of  the 

trajectory  to  the  horison.  As  already  aeationed,  flight  range  Li*  is 
not  included  in  the  auaber  of  paraaeters,  since  no  liai tat  ions  on  the 
flight  range  of  the  LV  are  inpoaed  by  the  conditions  of  the  probleas. 
The  faaily  of  disturbed  trajectories  will  have  a corresponding  faaily 
of  switching  points  Thus  we  can  plot  a certain  surface 

for  switchings  of  the  roll  angle  in  the  space  of  phase  coordinates 
V*  H,. , Generally  it  is  difficult  for  onboard  facilities  to 
"aeaorize"  such  4 surface,  and  in  order  to  use  it  we  aust  have  a 


DOC  » 0196 


PAGE  5 


computer  for  calculating  the  current  values  of  V,  9 , and  9 aboard  the 
LV. 


Let  us  use  the  very  widely  eaployed  aethod  of  transition  froa  v, 
9,  H to  certain,  to  sone  degree  equivalent,  parameters,  which  are 
aore  accessible  to  aeasureaent  by  siaple  onboard  facilities,  for 

t 

exaaple,  apparent  velocity  Vs  — go^  nxdt,  loadf actor  nx  and  the 
derivative  froa  the  loadfactor  with  respect  to  time  nx.  Obtain  lag  ns 
involves  certain  technical  difficulties,  and  thus  to  plot  the  control 
algorithn  we  aust  first  exaaine  only  two  paraneters  V^  and  nK.  la 
this  case  the  three-diaensional  switching  surface  Vint  H,n  breaks 
done  into  a certain  curve  in  the  plane  of  the  phase  coordinates  VSa 
aad  *,«•  Para  sets  r «,»  can  be  replaced  by  paraaeter  - switching 
tiae,  which  is  calculated  froa  a certain  fixed  value  of  loadfactor  nx 
■ n*  apparent  velocity  Vg  * 1*^. 

The  aain  factors  which  affect  selection  of  the  switching  point 
for  a LV  with  assigned  characteristics  aad  for  a certain  assigned 
entry  speed  are  the  scatter  in  altitude  of  the  conditional  pericenter 
of  entry  trajectories  *//„  iadeter Binancy  in  the  value  o t ataosphere 
paraneters  dp,  error  in  the  entry  velocity  of  the  IT  , scatter  in 
the  projected  ballistic  paraaetars  of  the  LV  APX  and  ^K,k,  , Me 
assuaed  aaxiaua  values  of  the  disturbances  which  are  givea  are  given 
below: 
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r, 


Navigation  error  with  respect  to  altitude  of  *+*dLt. ». 
per icen ter  $150  kw 


Error  in  valae  of  available  quality  AK.  .....  +100/o 


Error  in  value  of  reduced  load  on  frontal 


surface  APX +10°/a 

Error  in  entry  speed  ♦O. 2 kn/s. 


Oscillations  in  ataospheric  density  are  considered  within  the 
liaits  of  the  ainiaal  and  naxinal  aodels  [ 2 J: 

P«l.  = 0,012*-(u,w;  Pmm  = 0,013tf"oww,  Pm,,  = 0,01 9« ~0  07W. 

Za  the  case  of  a disturbance  of  any  oae  type  (for  exaaple, 
disturbance  with  respect  te  a//.)  we  sill  have  a certain  switching 
line  in  the  fora  of  the  depend—  oo  or  tor  — aeple,  in 

Fig.  1 for  a tf  with  Pf  » 80  kgf/a*  we  have  a fanily  of  switching 
lines  for  the  aaia  disturbances  (all  switching  disturbances  are 
assued  to  be  independent  with  a noraal  law  of  distribution) . Here  we 
see  that  the  position  of  the  switching  point  is  nainly  influenced  by 


fi 


fl 
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the  change  in  paraneter  la  ft«.  The  svitching  line  corresponding  to 
this  case  (we  call  it  noninal)  will  also  be  used  to  plot  the  control 
algoritha  (Figs.  2-4).  The  position  of  the  noninal  switching  line  in 
plane  n **•  Vsb  (or  *«.  Vj«)  nest  be  deter  wined  taking  into  consideration 

n 

M 

the  dynanics  of  notion  of  the  vehicle  about  the  center  of  the  nasses. 
Therefore,  after  we  have  selected  the  switching  line  on  the  basis  of 
analysing  noninal  trajectories  which  call  for  instantaneous  switching 
of  the  roll  angle,  we  nust  replot  this  liqe,  this  tine  considering 
the  turning  dynanics  with  respect  to  the  roll  angle  (see  Figs.  2-4). 
Exact  deternination  of  the  location  of  the  LV  within  the  navigation 

I 

corridor  of  the  entry  (Aft*  * eiOO-vISO  ka)  would  aake  it  possible  to 
further  define  the  switching  nonent  of  the  roll  angle  and  focus 
attention  on  other  disturbing  factors,  thus  assuring  that  we  obtain  a 
nininal  value  for  a final  velocity  at  the  assigned  final  altitude. 
However,  existing  net hods  for  independently  deter nining  entry 
conditions  (in  the  siaplest  forn)  do  not  allow  us  to  solve  this 

& 

problea.  Thus,  for  exaaple,  the  aethod  which  uses  inforaation  on  the 
nature  of  change  in  any  paraneter  aeasured  on  board  the  LV 
(loadfactor,  integral  fron  loadfactor,  etc.)  over  a certain  tine 
interval  [4],  produces  error  in  deterainiftg  ft.  on  the  order  of 


1 


150-170  ka  in  flight  near  the  lower  boundary  of  the  navigation 
corridor  of  entry.  For  this  reason  we  present  below  a synthesis  of 
control  without  specification  of  the  initial  entry  conditions. 
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The  studies,  part  of  the  results  of  which  are  shown  in  Figs. 

2-4#  indicated  that  the  noninal  switching  line  can  be  approxiaated 
with  sufficient  accuracy  by  sinple  linear  or  quadratic  dependences  of 
loadf actor  n**  or  tiae  on  apparent  velocity: 

^ (Cj  n 4"  bj  n -)-  aj)\  M) 


^ (,mJ  l^Jn  + kj  Kyn  + Ij), 

J 


' 


■her*  K» r.«  is  t 1m  current  value  of  apparent  velocity,  aeasured  on 
board  the  LV.  During  the  descent  process  the  current  aeasured 
parameters  <w  are  constantly  coapared  aith  the  calculated 

/S,CH-  As  soon  as  the  current  value  of  loadfactor  or  tiae  becoaes 
greater  than  the  calculated 


the  coaaand  is  issued  to  saitch  ^»<t>  froa  Kn>in(im.«)  to  /C„..x(t  = 0).  To 
increase  the  working  reliability  of  the  control  system  this  coaaand 
is  issued  with  a certain  delay  tiae,  during  which  the  conditions  of 
(5)  are  checked.  Then  there  begins  the  process  of  transition  of  the 
IT  to  flight  with  — the  duration  of  which  depends  on  the 

paraaeters  of  angular  stabilization  and  external  conditions. 
Thereafter  the  flight  of  the  LV  is  stabilized  with  until  the 

soft  landing  systea  is  activated. 

Evaluating  efficiency  of  control  algorithns.  The  efficiency  of 
the  control  algorithm  was  evaluated  as  follows.  The  main 
disturbances,  whose  aaxiaal  values  were  given  above,  are  assumed  to 
be  randoa,  independent  of  one  another,  and  having  a normal  lav  of 
distribution: 


* 
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//.=  //«  +4,  A//.  mi)  K = K + At\Kan, 

P x = P x + ^»  A/* ■•»>  v^„  = W',,  + Ai&Vn  m,I(  P=“P  + ^S^fnii- 

tor*  H.,K  9 PX9  ^IH  P rtpcunt  the  nominal  values  (aatheaatical 
expectations)  of  the  altitude  of  the  conditional  pericenter  of  the 
entry  trajectory,  the  quality  of  the  IT,  load  op  the  frontal  surface, 
entry  velocity,  and  density,  respectively; 

A//.  max*  AACmll,  A Px  max>  max*  Apm.x  — the  naxinal  value  of  the  expected 

deviation  of  the  corresponding  paraaeter  fron  the  noainal  value;  Aa_s 
- randoa  values  with  noainal  law  of  distr ibut iop,  zero  mathematical 
expectation,  and  aean  square  deviation  of  1/3.  Accuracy  in 
deteraining  the  final  speed  was  estiaated  using  the  aethod  proposed 
by  B.  G.  Dostupov  [5].  Here  the  value  of  the  aatheaatical  expectation 
of  the  final  velocity  and  the  triple  aean  square  deviation 

A Vt  = + 3s  were  determined  with  the  control  systea  working  and  with  the 
sub  effect  of  all  disturbances. 

Numerical  results.  The  effectiveness  of  the  work  of  the  control 
algoritha  was  estiaated  according  to  the  described  aethod  for  two 
types  of  LT  with  the  sane  value  of  available  quality  K,,tn  *0.3  but 
with  different  values  for  the  reduced  load  on  the  frontal  surface: 
Pj,= 80  kyf/a*  and  250  kgf/n*.  The  aaxiaal  value  of  the  rell  angle  was 
assuaed  eganl  te  t-«  = 180°  fee  IT  with  P*  = 80  kgf/e*  and  t««— 135°  for 
LT  with  Pm  — 250  kgf/B*. 


j 
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The  aoveaent  of  the  LV  (taking  into  account  turning  dynaaics 
•ith  respect  to  roll)  and  the  control  algoritha  were  calculated  on 
the  digital  coaputer.  Here  we  dealt  with  direct  entry  into  the 
Martian  ataosphere  (VIS“=6  ka/s) . The  ainiaal  peraissible  flight 
altitude  was  assuaed  equal  to  the  final  altitude 

The  effect  of  disturbances  on  the  value  of  the  final  velocity 
was  estiaated  for  two  variations  of  the  control  algorithm#  the  first 
of  which  uses  a switching  line  in  the  fora  of  dependence  nx.(Vs »)  (see 
Pigs.  3 and  4),  the  second  - in  the  fora  ef  tn(Vs«)  (see  Pig.  2). 

In  the  first  case  the  switching  line  is  approxinated  with 
sufficient  accuracy  by  a straight  line  for  the  LV  with  Px  — 80  k*«/e«: 

nx. Vs.  + 12,5 

and  a polynoaiaal  of  the  2nd  degree  for  a LV  with  P,  = 250  kgf/a*: 

In  the  second  case  the  switching  line  is  best  approxinated  by 
two  straight  lines  for  the  LV  with  P,  — SO  k«f/a>: 

wbea 


Vs<  1.M  ka/s  11.1  V„+.  a«.7f, 


r> 

i 
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•Am  Vs>  1.88  km/m  <.=  23.6  V,n+  0.2. 

The  values  of  the  final  speed  which  the  vehicle  attains  at  a 
certain  final  altitade  under  the  influence  of  different  types  of 
disturbing  factors  in  the  case  of  controlled  descent  using  the 
described  algorithm,  are  shown  in  the  table. 

The  natheaatical  expectation  M(VU)  end  triple  nean  square 
deviation  of  the  final  velocity  ai/k  — + 3a  for  the  coabined  action  of 
all  disturbances,  deternined  by  the  nethod  of  B.  G.  Dostupov,  are, 
respectively: 

for  algorithn  with  switching  line 

H.  - 3,5  km,  P?u  = 80  M{Vk)  = 246  y*/s,  A Vt  - 3o  - 

-47  rm/*; 

H,=>2  km,  P?H  = W M(VU)  = 208  *,/»,  A y,  - 3=>  = 

= 32  */•;  ' 

H,  = 5 km,  P?u  - 250  Aef /**.  M(V,)-  444  */ « A V,  — 3o  = 

= 42  w*/»;  ' 

for  algorithn  with  swtiching  11m  ntm(VSm) 

— 3,5  p;°“-80  VM  ^(^,)  = 237  *e/j,  A ^.=43  *./«; 

H.  = 2 KM,  P:OM  = 80  iff M (l/,)  = 206  *./5,  AV,  = 33  ««/*. 

If  we  coapnre  the  work  of  the  control  algorithn  under  "real" 


t 


* 
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conditions  (taking  into  consideration  the  turning  dynamics  of  the  L V, 
effects  of  various  disturbances)  with  "ideal"  control  (instantaneous 
turning  of  the  LV,  optiaal  switching  point  for  each  specific 
trajectory)*  then  we  note  that  possible  error  in  the  work  of  the 
algoritha  for  control  of  the  final  speed  is  (3«): 


on  « vehicle  with  Px  =^kgf/n*  and  Kv,c„  0.3 1 * 3.5  ka,  A V,  » 68 

a/s  when  tn(Vsn)\  H,  • 2 ka#  M/,  « 45  n/s  when  t„(Vso)',  ■ 3»5  ka,  A 

* 55  n/s  when  nxn(Vsn)\  H,  » 2 ka,  av,  * M a/s  when  "..(I'si); 


on  n vehicle  with  * 250  kgf/n*#  Kf,tn  ■ 0.3:  • 5 ha*)  Ai/»  ■ 

56  a/a  when  MV's*). 


Shown  for  coaparison  in  the  table  are  data  for  the  case  where  in 
Place  of  the  switching  line  a single  switching  point  is  used.  This 
corresponds  to  entry  along  the  center  of  the  corridor.  Apparently  in 
this  case  it  is  virtually  iapossible  to  achieve  reliable  descent  of 
the  LV  onto  the  surface  of  the  planet. 


These  aaterials  show  that  the  effectiveness  of  the  working  of  a 
LCS  using  the  proposed  algoritha  is  virtually  Independent  of  the 
aethod  of  asaigniag  the  switching  line  in  the  fora  of  or 

f.(V«a)-  The  eaiiaal  scatter  in  the  value  of  the  final  velocity  for  the 
working  control  systea  decreases  with  a decrease  in  the  final 
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altitude  Hx  «t  which  the  soft  landing  systea  Makes  its  entry. 
Methodological  error  constitutes  a value  on  the  order  of  50  a/s. 

The  proposed  control  algoritha  can  be  used  in  plotting  control 
systeas  of  final  descent  speed  for  vehicles  vith  an  average  value  of 
reduced  load  on  the  frontal  surface  and  vith  a low  value  of  available 
quality  at  various  rates  of  entry  into  the  Martian  ataosphere.  The 
algoritha  can  be  achieved  by  using  siaple  Measuring  devices  (for 
exaapLe,  the  acceleroaeter,  integrator,  clocks)  and  a siaple  coaputer 
for  arithaetic  and  logic  functions.  The  control  systea  can  be 
adjusted  by  a siaple  conversion  of  the  coefficients  which  approxiaate 
the  switching  line  and  are  stored  by  the  onboard  coaputer. 


Manuscript  received  2 January  1971 
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Fig.  2.  KB*:  (1)  s,  (2)  kgf/a*,  (3)  ks/s,  («)  Switching  Use  *m  (\) 
without  r dynamics  considered,  (5)  swtiching  line  ^6^,)  with  r 
dynanics  considered,  (6)  ks/s. 
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Pig.  3.  kBT:  (1)  Switching  line  without  y dynaaics  considered,  (2) 
ka*  (3)  Switching  line  with  y dynamics  considered,  (4)  ka/s. 

Pig.  4.  KEY:  (1)  Switching  line  without  y dynamics  considered, 
(2)  Switching  line  with  y dynawics  considered,  (3)  kgf/a*,  (4)  kn/s 
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Table.  KEY:  (1)  Value  of  final  velocity  under  influence  of  various 
disturbing  factors,  (2)  kgf/s*,  (3)  s/s,  (4)  ks/s,  (5)  Switching 
lines,  (6)  Switching  points. 
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